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Abstract 

We generalize the grid-projection method for the construction of quasiperiodic tihngs. 
A rather general fundamental domain of the associated higher dimensional lattice is used 
for the construction of the acceptance region. The arbitrariness of the fundamental domain 
allows for a choice which obeys all the symmetries of the lattice, which is important for the 
construction of tilings with a given non-trivial point group symmetry in Fourier space. As 
an illustration, the construction of a 2d quasiperiodic tiling with twelvefold orientational 
symmetry is described. 

Introduction 

The interest for non-periodic tilings first arose from problems in mathematical logics (Wang 
1965, Robinson 1971). However, since Penrose's invention of his well known non-periodic 
tilings (Penrose 1974, 1979, Gardner 1977), the motivation has changed to the study of geo- 
metrical questions related to such tilings. J. Conway (see M. Gardner 1977) and N. G. de Bruijn 
(1981) have played a dominant role in this field. 

We define a tiling as a covering of the plane by translations of a finite number of polygons 
with no holes or overlaps. Note however that, depending on the context, other definitions 
may be more appropriate, see e. g. Griinbaum & Shepard (1987). Here we are interested in 
quasiperiodic tilings. By this we mean: 

(i) The tiling is not periodic: There exist no translations (exept the identity) which leave 
the tiling unchanged. 
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(ii) If we put a (5-function to each vertex of the tihng, the Fourier transform of the resulting 
structure is a sum of (5-pcaks, whose positions are integer hnear combinations of a finite 
set of vectors {ki, . . . , k„}: 



(iii) Any finite part of the tihng appears infinitely often in the tiling. 

Condition (iii) is often dropped, but the tilings we will consider have this property. 

Quasiperiodic tilings may have symmetries in Fourier space which are incompatible with 
a periodic structure and do therefore not occur for crystals. A famous example is the Penrose 
tiling with fivefold symmetric Fourier transform. This has led physicist's interest to quasiperi- 
odic tilings, since icosahedrally symmetric diffraction patterns of an Al-Mn alloy, observed 
by Shechtman et al. (1984), could be explained in terms of a thrccdimcnsional version of the 
Penrose tiling (Mackay 1981, Duneau k. Katz 1985, Elser 1986, Kalugin, Kitaev k, Levitov 
1986, Levine &; Steinhardt 1984). More information on the symmetry of quasiperiodic tilings 
and the connection to physics can be found in the Les Houches proceedings (1986). 

The main concern of this paper will be a generalization of the grid-projection method 
used by various authors (Duneau k. Katz 1985, Elser 1986, Gahler & Rhyner 1986, Kalugin, 
Kitaev & Levitov 1986, Korepin 1986, Kramer & Ncri 1984, Socolar, Steinhardt & Levine 
1985). Our algorithm projects part of a n-dimensional lattice T in i?" onto an irrationally 
embedded d-dimensional subspace. It is based on a periodic tiling of by copies of a rather 
general fundamental domain of F, as opposed to Gahler k Rhyner (1986), who considered 
only fundamental parallelotopes. This extension allows for the choice of a fundamental domain 
whose closure obeys all the symmetries of the lattice, which is important for the construction 
of tilings with specified symmetry in Fourier space. The tilings so obtained are then no more 
tilings by parallelotopes alone, but can contain any kind of convex polytopcs. These tilings 
will not be the final step however. We rather prefer to consider their Voronoi partitionings, 
since these seem to have more relevance to physics (Jaric 1986, Henley 1986). 

The outline of this paper is as follows. After briefly reviewing the concepts of Voronoi 
partitioning and tilings by fundamental domains in sections 1 and 2, we present our generalized 
algorithm in section 3. In section 4 we prove that the tilings constructed in section 3 are indeed 
quasiperiodic in the sense explained above. Finally, section 5 is devoted to an example to 
illustrate these techniques. 

1 Tilings by parallelotopes and Voronoi domains 

The standard grid-projection method (Gahler k Rhyner 1986) yields tilings of the Euclidean 
space E"^ by a finite set of parallelotopes. Possible sets of paralellotopes can be obtained as 
follows: Let {ti}i=i,...,„ be a set of n vectors in {n > d). Any subset of d linearly inde- 
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pendent vectors of this set spans a parallelotope (see next section). With the parallelotopes 
obtained in this way i?^ can be tiled both periodically and quasiperiodically. 

It should be noted that the set of vertices of such a tiling is a special case of a De- 
launey (r,R)-system (Delauney 1937), which is a point set {vj} in E"^ with the following two 
properties: 

(i) The minimal distance between any two points of the system is r > 0. 

(ii) Inside or on the surface of any ball of radius R, no matter where we put its center, there 
is at least one point of the system. 

With each (r,R)-system we associate its Voronoi partitioning of E'^, which divides E'^ into 
Voronoi domains, also called Dirichlet or Wigner-Seitz cells. The Voronoi domain associated 
with Vj is a convex polytope and consists of all points of E*^ whose distance to Vj is not 
larger than the distance to any other point Vi ^ vj in the system. Note that according to 
this definition a Voronoi domain is a closed set, which means that boundary points belong 
to two or more Voronoi domains. An (r,R)-system has in general infinitely many different 
Voronoi domains, but a quasicrystal can have only finitely many (up to shifts), see section 3. 
A (periodic) lattice has even only one type of Voronoi domain. 

The Voronoi domain around a point Vj of an (r,R)-system in E'^ can be constructed as 
follows. Consider the set of all vertices Vi inside a closed ball of radius 2R centered at Vj. 
For each point Vi in this set, construct the {d — l)-dimensional hyperplane perpendicular to 
the segment Vj — Vi and passing through its midpoint \vi + ^vj. Each of these hyperplanes 
cuts into two halfspaces. The Voronoi domain Vj is the intersection of all those halfspaces 
which contain Vj. 

2 Fundamental domains of a lattice 

Consider a lattice F in E'^ generated by n linearly independent vectors ei, . . . ,e„. A funda- 
mental parallelotope Fp of F is the set 



Since a lattice has infinitely many lattice basis, it has also infinitely many fundamental par- 
allelotopes. A fundamental parallelotope is a special case of a fundamental domain F of a 
lattice, which is a measurable set with the following two properties: 

(i) The translates of F by all lattice vectors of F cover E"' with multiplicity one. 

(ii) F contains exactly one point of the lattice. 

It should be noted that a fundamental domain is neither closed nor open, only a part of the 
boundary belongs to F. In the following we will restrict ourselves to fundamental domains 
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whose closure is a convex polytope. Particularly interesting is a fundamental domain whose 
closure is the Voronoi domain. From property (ii) it follows that with each fundamental 
domain F a unique lattice point ^{F) is associated. This will become important later. 

3 The generalized grid-projection method 

Let us decompose the space E'^ containing the lattice V into two orthogonal subspaces, E'^ = 
© E-^. We will assume in the following that this decomposition is irrational, i. e. neither 
nor E-^ contain any lattice vectors of F. Let Fo be a fundamental domain of F, and Vr{Fo) 
the partitioning of E'^ into all F-translates of Fq. The closures of the F-translates of Fq will 
be called the cells of the partitioning. We assume that the cells are convex polytopes, and 
that the partitioning Pr(-^o) is face-to-face, i. e. if two cells have a non-zero intersection, then 
this intersection is a (common) face of these two cells. Note that the Voronoi partitioning 
is always face-to-face. If a cell is not equal to a Voronoi domain of the lattice, we moreover 
assume that the partitoning is generic in the sense that a face of dimension m is contained 
in exactly n — m + 1 cells. The dual of such a partitioning is simplicial, i.e. all cells (and 
their faces) are simplices. Next, consider an d-dimensional (afHne) subspace E of E'^ which 
is parallel to E^K We assume that E is located at a generic position, so that only faces of 
dimensions n-d to n-l of the cells have non-zero intersection with E. 

The generalized projection method now is described as follows. The vertices of the tiling 
are obtained by projecting orthogonally onto E the set W of lattice points whose associated 
fundamental domain has a non-zero intersection with E: 

W = {j{F)\FnE^ili,FeVriFo)}. (3.1) 

Next, we have to divide E into tiles by specifying all their faces of dimensions up to d — L The 
1-dimensional "faces" are obtained by connecting all those vertices by a straight line whose 
associated cells share a common face of dimension n — 1 which cuts E. If c? = 2, the tiling is 
then completely specified. For d > 2 however, the situation is somewhat more complicated. 
Those lattice points whose associated cells share a face of dimension k are the corners of a 
convex polytope of dimension n — k. For a generic partitioning this is evident, for there are 
always exactly n — k + 1 such points. For the Voronoi partitioning, we can argue differently. 
The points whose cells share the fc-face under consideration can all be connected by a chain 
or net of straight lines each of which is perpendicular to an {n — l)-face containing the A;- face 
and thus perpendicular to the A;-face itself. Therefore, all these points are contained in a 
single plane of dimension n — k perpendicular to the A;-face. Hence, in both the Voronoi and 
the generic case we can build the (n — fc)-dimcnsional polytope dual to a given A;-facc. If now 
the fc-face cuts E, we project its dual polytope to E. In this way we obtain a prescription 
for the subdivision of E into tiles. Note that with each projected dual of a A;-face also all 
its boundaries are projected, since the A;-face is contained in the corresponding {k + l)-faces 
which cut E too. 
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The same tiling can also be obtained as the dual of a grid G. This grid is given by the 
intersection of the Tinion of the boundaries of all cells of the partitioning with the subspace 
E. The grid divides E into convex polyhedral cells, called meshes^ the faces of which are the 
intersections oi E with the (n — l)-dimensional faces of the cells of the partitioning. Each mesh 
of the grid corresponds to a cell which cuts E. Therefore, with each mesh we can associate the 
projection of the corresponding lattice point, and two lattice points belonging to meshes with 
a common [d — l)-face have to be connected by the projection of the corresponding lattice 
vector connecting the two lattice points. The vertices associated with the meshes sharing a 
common A;-face will become the corners of a (c/ — A;)-face of a tile (these vertices are indeed 
contained in a [d — A;)-plane as explained in the previous paragraph) . In this way we see that 
the tiling obtained previously by projection can be reconstructed from the grid. According 
to this construction, it is the dual graph of the grid. 

What is not immediately clear is whether there will be overlapping tiles, i. e. whether the 
tiling is folded. Whether there are additional conditions required to avoid overlapping, and 
what these conditions would be, we leave as an open problem. For the Voronoi case however 
we have some (numerical) evidence that overlapping does not occur, and we conjecture that 
this is generally true for the Voronoi case. For the classical grid method, the necessary and 
sufficient non-overlapping conditions have been determined (Gahler &; Rhyner 1986, de Bruijn 
1986). 

From the grid picture and from the periodicity of F it follows that the tiling consists 
only of a finite number of different tiles (up to translation), for there are only finitely many 
inequivalent (n — d)-faces of the cells of the partitioning which can cut E (note that the type 
of such an (n — d)-face determines which vertices belong to the associated cell) . By a similar 
reasoning one finds that there are only finitely many arrangements of cells which share a 
common vertex, so that the Voronoi partitioning of the tiling, as constructed in section 1, 
consists of a finite number of different cells too. Bounds on the number of different patches 
of radius R of such a quasiperiodic tiling have been obtained by Gahler (1986). This number 
is finite and can grow only with a fixed power of R. 

Let us compare our construction briefly with the algorithm proposed by Gahler &; Rhyner 
(1986). They consider only special fundamental domains, namely parallelotopes. This has 
the disadvantage that for non-orthogonal lattices it is impossible to choose a parallelotope 
which is invariant under the whole point group of the lattice F. The choice of a symmetric 
fundamental domain is essential for the construction of quasiperiodic tilings which have the 
corresponding symmetry in Fourier space. By allowing a more general fundamental domain, 
e. g. the Voronoi domain, this deficiency is removed. This additional freedom is the main 
difference as compared to Gahler & Rhyner (1986). Using different "grid-" and "tiling-spaces" 
or including a subsequent linear transformation applied to the tiling could of course also be 
incorporated into the present algorithm. 
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4 Proof of quasiperiodicity 



In this section we demonstrate that the tihngs constructed in the last section satisfy the three 
conditions for quasiperiodicity formulated in the introduction. Since the proof of condition 
(ii) is a standard one (see e. g. Gahler &; Rhyner 1986, Zia &; Dallas 1985), we restrict ourselves 
to conditions (i) and (iii). 

First we prove non-periodicity. Let us define the projectors P" and P-^ projecting orthog- 
onally onto and E-^ respectively. Further, define the strip S as 



where the acceptance region M is the projection P^F onto with F a translate of Fq 
centered at E. Then, we can write the set W defined in (3.1) as W = T S. Clearly, 
the projection of W onto VF" = P^W, is the set of vertices of the tiling. Due to the 
irrationality of the embedding of and E-^, the sets P"r and P-^T are dense in and E-^, 
and there is a one-to-one correspondence between T, Pl'F and P^T, as well as between W, VF" 
and W-^ = P^W. Suppose that W^" is periodic, i. e. W^" is invariant under a translation 7||. 
Then 7|| maps vertices to vertices and is therefore the projection of a lattice vector 7. Hence, 
W^" is as well the projection of the set W + j. Since there is a one-to-one correspondence 
between W and 1^" this means that W = W + j. Let us project this equation to E-^: 
W'^ = W'^ + P^'j. Since the closure of W'^ is compact, this means that P-^"f = or 7 € i?", 
which contradicts our assumption of an irrational embedding of EK Therefore the tiling is 
non-periodic. 

Now we show that every finite part Wj C W^" has infinitely many copies in W^K Denote 
by Wf the unique subset of W such that Wj = and by its projection to E-^. 

Since is at a generic position, no lattice points are projected onto the boundary of M, and 
so Wf - is in the interior of M. Let A be the distance of Wf - to the boundary of M, 



For every lattice vector 7 whose projection onto E is inside an open ball of radius A, we 
have that the finite set Wf = Wf + 7 C F projects into M, P^Wf C M, and therefore 
belongs to the strip S. This means that a translation by P"7 maps Wj onto an equivalent 
set. Since P-'-F is dense in E-^, there are infinitely many such lattice vectors, and so the proof 
is completed. 

5 Example: a dodecagonal tiling 

As an application, we discuss the construction of a class of twodimensional tilings with twelve- 
fold symmetric Fourier spectrum. These tilings have first been constructed by Stampfli (1986) 



S' = {m-|-e|mGM,eG E^, 



(4.1) 



A = min (|x — y|). 
y&dM 



(4.2) 
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by means of a grid. They might be relevant for the description of quasicrystalline Ni-Cr (Ishi- 
masa, Nissen & Fukano 1985, Gahler 1987). More details about these and related tilings can 
be found in Gahler (1987). 

The relevant lattice for our case is the diisohexagonal orthogonal primitive lattice (Brown 
et al. 1978) in four dimensions, denoted by T. This lattice has a point symmetry group which 
contains the subgroup D24. The latter is the relevant symmetry group for our purposes. The 
lattice r is easily constructed as follows. Let us decompose into two orthogonal subspaces, 
E'^ = © E-^. is the space onto which we will project. Let {ei, . . . ,612} be a star of 
twelve vectors in E^ such that their projections onto E'" and E-^ are given by 

e'l = (cos(7r(i - l)/6), sin(7r(i - l)/6)) 
e^^ = (cos(57r(i - l)/6), sin(57r(i - l)/6)) 

with respect to Euclidean coordinates in E'^' and E^. The vectors {e^, . . . ,612} can be ex- 
pressed as integer linear combinations of the remaining four vectors, and since {ei, . . . ,64} 
are rationally independent, the set {ei, . . . ,612} generates a 4-dimensional lattice which will 
be identified with F. 

In this example, we choose a fundamental domain whose closure is the Voronoi domain. 
Therefore we have to construct the Voronoi partitioning of T. We note that the space spanned 
by ei and 63, denoted by E"', is orthogonal to the space E^ spanned by 62 and 64. The 
vectors ei and 63 generate a 2d regular hexagonal lattice F'* in and 62 and 64 generate a 
corresponding lattice F^ in E^. Therefore, F is given as an orthogonal sum of two 2d regular 
hexagonal lattices, 

F = F"©F^ (5.2) 

Next we recall the fact that in such a case the Voronoi domain of F is given by the topological 
product of the two Voronoi domains of F" and F^, 

V = V'' X V\ (5.3) 

which of course are regular hexagons. Let and i?^ be the hexagon nets given by the 
boundaries of all Voronoi domains of F" and F^ respectively. Then the union of the boundaries 
of all Voronoi domains of F is given by 

iV = (i7» X E^) U (E" X H^), (5.4) 

i. e. N is the union of two orthogonal arrays of hexagonal "tubes". Let E he a generic 
plane parallel to EK The grid, i. e. the intersection of A'^ with E, is then the union of the 
intersections of the two arrays of tubes, which are both regular hexagon nets, turned with 
respect to each other by 30° (see Fig. 1). The elementary hexagons of these nets are by a 
factor of two larger than the projections of the Voronoi domains of the hexagonal lattices. 
The relative positions of the two nets are determined by the position of E. 

For the construction of Stampfli's tilings, the two algorithms discussed in section 3 now 
read as follows: 
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Fig. 1: Grid given by two hexagon nets 



Fig. 2: Quasiperiodic tiling dual to grid of Fig. 1 



A: Projection construction 

Project the center of all those Voronoi domains of F onto E^^ which cut E. Connect all 
those points by a straight line whose Voronoi domains have a face in common which 
cuts E. 

B: Grid construction 

With each mesh of the grid N OE, associate a vertex of the tiling. If the meshes of two 
vertices have a common face, these vertices are connected by a line of unit length which 
is perpendicular to this face. This is Stampfli's presciption. 

A tiling constructed in this way is shown in Fig. 2. 

Let us finally note that this is a particularly simple example, due to the fact that F is an 
orthogonal sum of two 2d lattices. From this it follows that the grid NCiE is the union of two 
simple periodic grids. In the general case, NDE would be very complicated, but nevertheless 
all our constructions would go through as well. 
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for the warm hospitality during the time when part of this work has been done, and to the 
USSR Academy of Sciences for financial support. F. G. would like to thank Peter Stampfli 
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